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mmm 

M^ffc cd ^ n -fe y -9- ^ a - ;b tc, *i ti * iT £ ^ n v V & taff-r -5 

i tcaa«©aH»«*. 

[»*^ 3 ] t& ffilKR © -f n -fe >y -9- =e a - ;v &m -r * /M?#I ^ n v V <d 
*y h&^x-T^y ^ 
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[f&#£6j tilted is W*&l&ni/x*-ftMlZ. MfB^J 

© 3^ e> S» W * r u x a ic «fc o t JE&ff #f © J? «f & ff v ^ , M s b*t #i © 

{C J: U , mm±& tC o ^ T L U £ T f 5 £ & #^ £ -f S » 2 tc ffi 
[0 0 0 1] 

o 

[0 0 0 2] 
[^*©&ffi] 

1 & £ ff J¥« tc J: o T 8? < JC t±, S£:£ 1 & #g5£ & L 
[0 0 0 3] 

. LU^» (LU Factorization) V^ ? #i£iC <fc fttf, « *±Eift 
?T#I (upper-triangular matrix ) £TH.^ffT#J (lower-triangular matrix ) 
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Cholesky Factorization) £ ^??T#f £ 0 

.^hMfl*^*fcfc^ny^S:LU^ (1) bfctl, #j£l,fcfT'<* 

y*&H*rLT (2) frfb, IE#/M??'J£3!§t (3) -TSfltJl&ifty 

[0 0 0 4] 

(l) ©^Ste, — oo^n-fey-y-CjS^JcfToTv^fc. 
^■t) SS©itiKi^£vvfB£:LTv^ 0 r©MJH (2) (3) tf)g|5#(Z>3g 
[0 0 0 5] 

m*J^Xh(D*£^ (3) <Z>fH*T*, <B#1 2nmt.'bz <X%. *«CZ>i^ 
#&#&ofc 0 M^^'JS^A51M*« (SMP) flltf/J^V^'ft 

[0 0 0 6] 
[0 0 0 7] 

b ) jew^t^©^ i fcjjm&amm 
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cyclic ic#^nizy-!Hc#ifcUTiMSS'tt, #7n iz v 

[0 0 0 8] 

[0 0 0 9] 
[0 0 10] 

[M&#8Si-rsfc«)©#ia] 

[0011] 
[0 0 12] 
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o 

[0013] 

&m* z evmxju?&#m$f-mi*. w&v-iumwi o - 1, 10-2, • 

••10 - n#2&3f * y S/a^y 13-1, 13-2, •••13-n 
LT^SS^^ 1 2 lC«8ffi3*l«. 1 0 - 1 , 10-2, ---l 

0-nli, •^■tf)ftg|S&& V^£, 2 V V V =l * *E V 1 3 - 1 , 13-2, • • 
• 1 3 -n <fc U^nizy-IHIUC l#C=3f^y S/rL*^ U #tgtt£>*l<5 0 #y 
nt»;tlO-l, 10-2, • • • 1 O-nlcM-tftotV^^tU^^a 
11-2, • • • 11-ntt, i5^il2^bT^nt»y 
-9-10-1, 10-2, • • • l o - n^T^irXBT|g^^oT^) 0 :/n-feytJ- 

io-i, 10-2, ••• i o - n^-^mm^nom^icit, 

U ^S^rL-^l 1 - 1, 11-2, • • • 1 1 - nfr£> 1 ■D<D7U J bv-V--tfU : % 
-t&T-#$:2?k*Vy*s=L*^V 1 3 - 1, 13-2, • • • 1 3 - n IZ&ffi 
U JglC, 2<^^-V^S/rL^^U^^^ ! a#fili:^-5T r -^ l^^r^J/a^ 

^y izn^-Lxmm^noo 

[0014] 

^*^i»S4x, 2t^^>yi/a^=E'Jft0f-^MTMIb»5i:, ^^'J 
11-1, 11-2, • • • 11 - n(D|*K mmz?- * £f#oT£ 

m&izit. _kM L £ J: e> (c, ;*^y ^^-;i/j^#^nizyi^S^£#^ 
-5 s - ^&2i^^t'yi/i^ : E'J iC|&M U l^dp-wi/a^^y ic^a-^©^ 

=ey ^^j.-;wc#£&», &©4M<z>fc#>jc, Stf^^u^^a-^^f- 
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fli(D-J'n-\zyV'<DW!im<DMffl : &:£Z>Z. £ S^'J Tl^il (Barrier S ynchronizati 
on) V>5 0 
[0 0 15] 

L 1 ~L 3 ©BP^tt, :/n y ?ipi£/h£ < bt, lo©^nt-yt (P 
[0 0 16] 
[0 0 17] 

7 s - WIS &£HR PEtH< # P E -e-r-*- * £J5ffift bti 

-5 z: £ >«M 4: * £ . 

[0 0 18] 
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[0019] 

X (b + k) X 8~8Mb y t e &ifet%©&^ny ^jSilbtfilfflt*. 
[0 0 2 0] 

[0 0 2 1] 
[0 0 2 2] 

o 

*^^n v^fiSr^i-s. #cic, *^y 2fC;fev^T, #;* i^y 
F^lt*^Q»;i'Sft)IU ^l/»;Ktffilt*^D«;^ (HI 2 ©D&tf 
Li) Srf^|S«{C3t:-'tS. fit, *7^y -fS 3 iCfc^T, #Xl^yKT*b? 
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[0 0 2 3] 

^Sfc/WT-f ;|&fc>*JT*&v^-£&Ctt, Xf^^S 5lC;l3^T, #XWyK 
"Cdfe^JtC^D >yHL (E|2#flg) ^ot, ^D^Ui (02« ^LU 

y KT'MJlC yn^y^Ci (SI2#Rg) & "/a v V L i i: -fU y V U (Ui5: 
mfriafr^t&O) : g2#I) <D«T*3g3fr"t£ 0 fit, J^y^S 2 
?^®yn>^ (C i ftM^S^n-/*) |S|^(D^f^T'LU^-rS 0 fit 

, #^fc/jN$ < fts^Ji^n y ? ^*^(c^n >y ? d &c;#jSi-£g|5#tf)#K:& 

[0 0 2 4] 

o 

^ZlXte, 2048X204 8tf)?f#J£4 ^ Wy FT L U##?-f £^-£-£#KC 
[0 0 2 5] 

£1% 2048X204 8 (DftFiZ^U V #(lg2 5 6 L U^^f -5 t><£> £ ~T 
£o y PCDmfrmt. ® 5{C^1-££*JT*&£o l^T, 4o©CPU ( 

Xl/>y F) •e«!aS:llfT , rS»^, FT?3*KfciK« ( (2 5 6 +4 4 8 

) X 2 5 6 : 8MB (L2^ft-yJ/a©^tS) J: »J/hS SrflS«U 0 

6 IC^-TJ: c) iC, #X l/y K0>#fB«tCD 1 + L 1 . D1+L2, D1+L3, 
D 1 + L 4 $:3tf--rs o 
[0 0 2 6] 
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ADH L/T, 



[0 0 2 7] 



y(n 3 M THREAD) x y~ 



= nb 



[0 0 2 8] 

nb^5124^ yn^^ii(g=5 12 
nb^2 5 64?), y #iflg= 2 5 6 

nb^l28^f), ^0^^(1=128 

fn«^, :/ny 6 4 

[0 0 2 9] 

lT&fc>*>, LU»n^h(t 2n 3 /3 {n\tmWkm T*, n 3 iCjt 
ffilfS. i#©37M:#THREAD"elJ!|ftU ■£© 1 %#* 

[0 0 3 0] 

7iC^-f 0 EJ7{C£V^&, LT = D1+L1+L2 + L3 + L4©§P^I:LU 
^t57;i/=f'JXA$:^bT^5. LTtt, 2 0 4 8 X 2 5 6 y 9 £ & 

[0 0 3 1 ] 

£t\ (1) ©BP4HC;J5V>T, bfilCy h Z$k%.-?Z>o iblks liLT©ipg-C*»J 
, 2 5 6-e&3 0 lengli, LTC^STfifey, 4*<£>^lir\ 2 0 

4 8T*&£ 0 j j Kite, tf/Ky hfiD^-r TMPlCtt, KT/Jfy h<D 

[0 0 3 2] 
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*OT, (2) ©«P^ICfeV^T, Slfc&SLTV^ LTft©?!I## i «fc U> tf 
t 1 -* ilAtlSA-Sc (3) <Z>g|5#T% #|i ©tT)J?v h ^lot LU^ 

[0 0 3 3] 

Z.(D (1) ~ (3) £i#l~iblks ttfcoti'JIliltS. 

LTCDI^T'&S, leng#=feo Zl *l £ ©MStt L 2 

vr=.$.£mm*fit> #*i^y Rc£tt3 7;i^';xA&0 8fc^-r 

[0 0 3 4] 

08&C;i3^TkL LTU±, n-#;Hgc, pivot(4), GPIVOT, ROW (ib 

iks ) it, ^ffl^ctci&Msns^-^-e&s. 

£1% (4) tcfcv^T, #Xl/«yFt'fcfjK«yh^5 0 fit, (5) T*** 
ffi£IB#Jpivot(4)©x Wy K#^©Mifc^t-5 0 (5) ©fttc, AUTRJB 
&^oT, (6) IC&V^, **fcT3j?y h Wy K##&GPIVOTlC|&|ft"t 

s„ fu, (6) <d&x% ntf^vrmmitnLz* (7) k^t, * 

^UrWJeSrflJCSo (8) &C£VMTk]:, GPIVOT^ 0©££tt> *^^*;bttD 

1 ©rtic»s*\ A*i#*.>Fin?&y, n-#;i/«K:Aft#a.£. gpivot##t 

H R E A D IC^£ L il -*&t>?>* GPIVOTtf 0 «fc »; V^-^JC tt, Jft^CbTjK 
yhS:j$fc355 ^^Uy Kli, R0W<£>p*3^£ i frBCDff^^ h;i/£A*l#;L£ = -£ 
IT, (9) , (10) IC£V*T, LU»feJ!)(DSI^^ fit, ±IH 

(4) ~ do) zmm-?z>-7nv ?<D-£x<Dmz-D^xfto 0 ^^t>*>, i~ 

iblks $.X(D i iZ^^^Tj&M&ifiikVM'fo 
[0 0 3 5] 

lit*, LT i (OLU^fttt'^U TlTOfciKSo #Xl/«yK 

®Dl®MttLLi:UUlCLU^»aftt^S. fit, Vv Kt'Ul*- 
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LL _1 U1, U2«-LL -1 U2, U3<-LL _1 U3, U 4 «- L L _1 U 4 1/ 
y KT*»JlCftSi:-r5 0 ZlOffW^^T*, Dl, LI, L2, L3, L4SO- 
# A IC 3 fcT- A y * U Ay rira^m^,, MIC, C1«-C1-L 
1XU, C2^C2-L2XU, C3«-C3-L3XU, C4«-C4-L4XU 

[0 0 3 6] 

E9jCa^StiS J:^tc, iiiB&JS & "T -5 £ i: tc £ o T , f?*J©*M(ff(Z>f?;&tf?!J 

-rj:e>{C, Dl, LI, L2,. L3, L 4 © -?U y V lC#tfT, #^l/^KJC3 
2 5 6 X 2 5 6®^ny^jSS. iCDfflJ^tt 1 o©^ 1/ y FtLU^jJf btffl 
[0 0 3 7] 

_tIH^-e«, LTiSLU^HtSiit, >>:3.*=EyJttf>7 i - 
# <#jM^££&i?!j§ift&Lu#^£^^^ ci^ci 

-Li XUCD^Jfli, ^^^^a^^U^^^-^&^tCflJfflbfe^ras^SE^I 
[0 0 3 8] 

Ell 0&, S»WLU^jg7^J'JXAS:Siit«0Tf!fc4. 

«t<&3E»tt, @9©r;i/JuxAT'ffitt^ lti (#xi/7FtDi + Li 

, k (LTi©1^7Ci©*S3) , iblks (^nyi'fi) ©fffiiC, L 
[0 0 3 9] 

ii-r, -9-^>-^>©«IOT, nwid#8, "f L U#»&*T^*S# 8 T? 

feSs^^^&fllWT^So YES -nlCfi, LTi (ist : k, ist : ist +nw 
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id-1) £M^CLU##?lf-£= ist : k^W^^gBtt, IBRtfist 

k £"e©L T i ZTjk^MVfcV, ist : ist + nwid- 1 £ V 1 * 3 ^Ifc^ist 
;fr£>ist +nwid- 1 £T*<Z)L T i *5fcT?*« 0 JBJtTlCfc^Tfcl^BSTf&S 

a 

[0 0 4 0] 

LT i ©LU^»fCfeVATli, _L3S (4) ~ (10) 0>#l3lS:f=f 3. fc£U 
?ftf>A*l#X.3P£M±, «£iblks T% LTi (i. 1 : iblks ) £A*l#;l£o 

_hffi2ffWT3^NO©»^tCt±, LU 3 *Bnifid& 2 T?fflo fete© L 
U^iK&^LU^^^^-^^SrSlliWtCl^tfffiLTfT^. TRS 
fcV^;i/~f- ^<Z>;i/-^>&, LTi (ist : ist +nwid/2 - 
1, ist + nwid/2 : ist +nwid) ^Jgif^-So HfC, LTi (ist : ist +nw 
id/2-1, ist : ist +nwid/2 - 1 ) (DT^ftftfUL L fcfljjg It, LL" 1 

LTi (ist H-nwid/2 : k, ist + nwid/2 : ist +nwid) =LT i (ist + 
nwid/2 : k, ist + nwid/2 : ist +nwid) —LTi (ist +nwid/2 : k, 
ist : ist +nwid/ 2 - 1 ) XLTi (ist : ist +nwid/ 2 — 1 , ist+nwid 
/2 : ist + nwid) 

b ) lEmttibftm&itkft i &#e^tf>«?& 

[0 0 4 1 ] 

Rc*f£[fiW £3g#rt-£#J^a y *9S#L l , L2, l 3 Sr^|g«jc 
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6. C2£C5, C3 £C4 *m&&t>l2:Tmm*:fio &olZ.-tZ>) . 
Ktti#B2fctf 2 X#T + 1 - i&E<D-fUy9&-X.ffi"r&Z.£T*&ffi&mmz. 

[0 0 4 2] 

IR&^XT V * S 4f #J#IE«*f fMrfllT? & 5 CD T*, i2©U Jc*flS -T -5 3P# tt, 
Ll+L2 + L3^^j;5^D«yi'0SlL T & o T tf>T% Z<DWf£lZ 
t±, L T fflS^tt, SlirSrff e>^M#M<, LI, L2, L 3 &f *if *ifg«lT 

[0 0 4 3] 

£1\ (1) ff^J©— *&<Z)^ny^S:=lU^#-^J!¥L, 3 0 CD 

3 1 lC3fcf-*r<6o fit, 3 1 C0T©^-eH^4xTVN^gp^$:^S 
^&MV>T5Sif-r^o (2) (C^VNT, ^St--5^JC0ifia5: 2>^tCLT 3 2 

CDgP#&3 3tC:3fc 0 -U 3 3 COT©^mT'H*tlfeg|5^$:MaBl5^&MV^T, 

fit, (3) lzm/uX% 3 4 3 5 (C=rt°- U 356DT© 

j^j»T-^snfeap^&, aa*&as#&^T35*H-s. xk, (4) ic^^t, 3 

6T'^$ti^§P^^3 7 ten if- I, 3 7 ©T©^" fc^THSr-r 
JSIC, (5) {C&^T, 3 8 £ 3 9 tcntf- I, 3 9 ©T©,£J&"e*t£ ftT 
v^gp#£^3P#&Mv^H§TU (6) &c£v>T, 40^:4 ncnif-l, 
*MK8P#&JBV*"C, 4 1 ©TCaP^'SrHSf'tS. (7) JCfcV^, 4 2 CD 

ap#£4 3lCnfcT-U #ttteP#&MV^4 3©T©S#&3€*r'*<&. ^©£3 

[0 0 4 4] 
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[0 0 4 5] 

£-f, &7l/y 4 IC^ftSDx&tfL i &3ii&fl&fc1!B«LT i fC 3 

tT-fS. fit, LTi^LDL^in. LDL T #8?&#!llf lft&#&t 

ffe>. fit, d l i ^&>TCDniw.vmzmmzzi 

[0 0 4 6] 

DLi<-DXLi T , dZLT% DttDx©^Mlt'fe l J, 

•fit, Dx H©L i ttM^Jtc, &£(Z)«c&c:ilf-A 

*; i'tfi. fit, ^t, ^UrraSB&5t»J, EI 1 5 ic^3*i£ J: o fc, Cl 
fcC8, C2£C7, C3hC6, C4 hC 5 &^7tlt, #XbyFtM?lJ 

•^l/^Flt tt, 
C 1 «-C 1-L1 1XDL1 1 
C 8 «-C 8-L42XDL42 

• 7 Vv F2 tft, 

C 2<-C 2-L12XDL12 
C 7 <-C 7-L41XDL41 

• * 1/ y F 3 t tt, 

C 3 «-C 3-L21XDL2 1 
C 6 «-C 6-L32XDL32 

• 7Vy F 4 t tt, 

C 4 «-C 4-L22XDL22 
C 5 «-C 5-L31XDL3 1 

[0 0 4 7] 

fit, r^*-e©sii¥*««F*)ofc^-e^uri3jes:«t*j, lJS»J/J>3<Ja: 
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FT*, LDL T ftMLTMmzmT-?Z> 0 
[0 0 4 8] 

V \*WL%: #THREAD£l/T\ L$:2X#THREADl{:^|!|b, C©TH 
ftBft*imV< 2 X #THRE ADflJC^-H-r*. fit, ±fcT*6#THR 

So 

[0 0 4 9] 

H16H:, LDL T ftM<Dnm&ITJl3 V XAS^tBt**. 
LDL T ^i07;^'JXAH t^-f>LDLi:lt||3Sa*l6. if?" 

[0 0 5 0] 

£t\ nwid^8«^lCtt, (2 0) <B«S#T?, tt® L D L T #8? fcfr 5 . ^ 
tt, L T i (ist + 8 : k, ist : ist + 7) SKSt"*"*. Z(Dt^. LTi 
(ist : ist +7, ist : ist +7) CD_hH^|g|J^}C D L T ^Aot^5©T% 
(DL T ) £IC ioTJIfft^n 

[0 0 5 1] 

(20) fflSUfflIFit% nwidtfS^J&V^^WTSflfcl&^Cfci:, nwid/2 
&§T&&nwid£: Utt^-f > L D L &Pg*tffflU ^T^S* HI', LTi 
(ist : ist +nwid/ 2 - 1 , ist +nwid/ 2 : ist +nwid- 1 ) {C D L T £ n 
fcf — 16„ DJ£, LTi (ist : ist +nwid/ 2 - 1 , ist :ist+nwid/2- 
1) <Dttfi&M"e%>*)* L{£, LTi (ist +nwid/2 : ist +nwid- 1, ist 
: ist + nwid/2 - 1 ) T*&U> i©L^iIt5 0 
[0 0 5 2] 

fbt, LTi (ist + nwid/ 2 : k, ist + nwid/ 2 : ist + nwid- 1 ) £ 

LTi (ist + nwid/ 2 : k % ist + nwid/ 2 : ist +nwid- 1 ) = L T i (is 
t +nwid/2 : k, ist +nwid/2 : ist +nwid- 1 ) —LTi (ist +nwid/ 
2 : k, ist : ist +nwid— 1 ) X L T i (ist : ist +nwid/ 2 — 1 , ist + 
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nwid/ 2 : ist +nwid— 1 ) 
"5. 

[0 0 5 3] 

LD L T »t^-f>LD L £?l!MJClWffi-r o fit, MS 
*5fiWfi!)m^»or^J 'J XAIi, CD-ROM^©Br«SB»8SWM>, ^-K 

^ a & y n -fe y -9- jc a - k l x m m -f z> mm ■& # tl h *i s „ 

[0 0 5 4] 

• LU^i©ittt, 1) > 2) , mmnUX*r-frM(DJ:Ml l t2) 

1) P.AMESTOY, M.DAYDE, and I. DUFF, "Use of computational kernels 
in the solution of full and sparse linear equations" , M.COSNA 

RD, Y.ROBERT, Q.QUINTON, and M.RAYNAL, PARALLEL&D I STR I BUTED ALGORI 
THMS, North-Holland, 1989, pp, 13-19 

2) G.H.Golub, C.F.van Loan, "Matrix Computations" , second editi 
on, The Johns Hopkins University Press, 1989 
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[0 0 5 5] 

m tr y v <d & e» mmm <d n -t y v ^ *? a. - n, ic &m. <d * ^ u « tc as 

[0 0 5 6] 
[0 0 5 7] 

T* & 5 ^ £ £ £ i" £ fB lie |B« © IB^^flSo 
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[0 0 5 8] 
[0 0 5 9] 

[0 0 6 0] 
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[0 0 6 1] 

*. 

[01] 

[02] 
[03] 

*5§^©H»ffiiC^o L U^»©M3?!I«ia©«E^S:SiWrSH 2 ) 

[04] 
[05] 
[06] 
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iblks 



leng 




DO i=1. Iblks 

TMP=0.0 DO;jj=0 
DO j=i. leng 

IFCABS LTQ. I)). GT . TMPJTHEN 

TMP=ABS(LT(j, 0) 

ENDIF 
ENDDO 



IF(jj. QT, i) THEN 
DO k— 1, iblks 
TMPX=LT(i. k) 
LT<i, k)=LT(jj. k) 
LX(jj. k)=TMPX 
ENDDO 

END IF 



(1) 



(2) 



DO k=i+1. iblks 

LT(i. k)=LT(i, k)/LT(i. i) 
ENDDO 

DO k=i+1. iblks 
DO l=i+1, lens 

LT(I, k) = LT(l, k)-LT(l. i) x LT(i, k) 
ENDDO 
ENDDO 
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ENDDO 
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ibika 




LTI 



DO i=1, iblks 

TMP— 0>0 DO;jj=0 
DO i=l . t«ngi 

IF<AB3 LTi(j. i)). GT . TMP)THEN 

TMP=ABS(LTi(j. 0) 

EN OIF 
ENDDO 

pivpolC #THREAD) =ii 
(#THREAOIi. *U^K0##, 4XU<yK 



] 



(4) 



(S) 



BARRIERWM 

iF(#THREAD. EQ, 1> 

jx=0:QPIVOT=0 
DO bc=1. 4 

IFtpivot(ix). GT. ix. AND. PIVOT(ix). GT. iblka)GPTVOT=ix 

ENDDO 
END IF 
BARRIER F5IIM 

IF(**THREAD. EO, GPIVOT)THEN 

IFCii, GT r Othen 

□O ix=n, ibBca 

ROW(ix)=LTi(ii. ix) 
ENDDO 
END IF 

BARRIER Plffl 

IFCGPIVOT. EO. OTHEN 
IF(jj. GT, DTTHEN 
DO 1 = 1, iblks, 
TMPW=LTi(i. ix) 
LTi(i, ix)=sLTiCii. ix) 
LTiCjj, ix)-TMPW 
ENDDO 
END IF 
ELSE 

1F(#THREAD, EO. G PIVOT) THEN 
DO lx=1, iblks 

LTiCii, ix) = LTi(L Ix) 
LTi(i, ix)=ROWCix) 
ENDDO 
ELSE 

DO tx=1, folks 
LT i(i. ix)=ROW(ix) 
ENDDO 
END1F 



DO k=l+1. Iblks.. 

LTiCi. k) = LTi(i. k)/LT(i. ») 
ENDDO 

DO k=i + l, iblks 
DO 1=1 + 1, lensi 

LTKl. k)=LTi(l, k)-LTia I) *LTlCl. k) 
ENDDO 
ENDDO 
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[01 0] 



subroutine LLKLTi, k, iblks, ist. nwid) 
(LTi{iaxU*yKl?D1 + Li*lMnV 
k : LTi<Z>T2K5t § fl)^*. 
iblks ;7ac£ ft 
ist:Lu$}«£&tf>«<&a. 
nwid:Lu»S£fr3rtl) 
IF(nwid, eq, ftlThan ( rtl8*<»/IO 
LTi(ist:k, ist. ist4-nwid-1)*^JI^LU»»-r« 0 
f" ZCt, (4!-(10)^ o 

fcf£LfT(OA*i**«B^(***ibik-riLTi(i. 1 , iblks)£A*i*;L£ 

else 

call LLKLTi, k, iblks, ist nwid/2) 
calf TRS< ) 

LTi(ist:ist+nwid/2 — 1, ist+nwid/2:ist+nwid) 
$M*rT-5oLTiCist:ist^nwidy'2-1 l ist: ist+nwid/2- 1 ) 

(DT=afr9JLL£*JJBLT* LL + *S^e,»ltrS«r 

call MM( ) 
f~ LTl(i8t+nwid/2:k. ist +nwid/2 : ist + nwid) 

a* LTi < ist+ nwid / 2 : k. ist + nwid/ 2 ; i st+ nwid ) 
— LTl(ist+nwid/2:k, ist: ist+nwid/2- 1 ) x 
v_ LTi(ist:ist+nwidy^2— 1, ist+nwid/2:iat+nwid) 

aamerfg»W 

call LUCLTi, k, iblks, ist+nwid/2, nwid/2 
end if 
return 

end subroutine 
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subroutine LTD(L_Ti, k, iblks, ist, nwid) 
IECnwid, EQ, 8)THEN (f|)8*<S/J*) 
DOi=ist, ist+7 ^\ 
DOj=i + 1. ist+7 
LTi(i, j)=LTi(j. i) 
LTi(j, i)=LTi(j ( i)/LTi(i, i) 
ENDDO (20) 
DO jy=i+1, ist+7 
DO Jx=Jx, ist+7 

LTJ(jx f jy)=LTi(jx. jy)-LTiGx. i> xLTi(i, jy) 
ENDDO j 
ENDDO J 



c 



LTi(LTi(ist+8;k. ist: ist+7) £5£«r 

LTi(LTi(ist:ist+7. 1st: ist+7) <D± = ft |CDL T ^A^TL^<DT? 
(PL T ) ^t>*MtTM»rt£ 



ELSE 

call LDL(LTi, k. iblks, ist, nwid/2) 

-LTifist ist+nwid^-r ist+nwid^2 : ist+nwid- 1> 
IZDL T $3tf— T^o (DttLTi(ist:ist+nwid/2-1, ist:ist+nwid/2-1 ) 

LTi (ist+nwid/2: ist+ nwid — 1 . ist: ist+nwid/2— 1 ) 



•LTi(ist+nwid/2:k ( ist+nwidy2: ist+nwid— D^Mft 

^ LTi(ist+nwid/2:k, ist+ nwid/2: ist + nwid— 1 ) ^ 
= LTi (ist : ist + nwidX2 : k, ist+nwid/2 : ist + nwid — 1 ) — 
LTi(ist+nwidX2:k, ist : ist + nwid— 1 ) x 

^ LTi(ist:ist+nwid/2— 1, ist+nwid/2 : ist+nwid — 1 ) 

CALL LDL (LTi, k, iblks, ist+nwid/2, nwid/2) 



ENQIE 
RETURN 
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